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D I Abstract. The ABC model is a driven diffusive exclusion model, composed of three 

species of particles that hop on a ring with local asymmetric rates. In the weak 
asymmetry limit, where the asymmetry vanishes with the length of the system, the 
2 ' model exhibits a phase transition between a homogenous state and a phase separated 

"^ ■ state. We derive the exact solution for the density profiles of the three species in the 

hydrodynamic limit for arbitrary average densities. The solution yields the complete 
phase diagram of the model and allows the study of the nature of the first order phase 
transition found for average densities that deviate significantly from the equal densities 
point. 



PACS numbers: 05.50. +q, G5.70.Ln and 64.60. Cn 



1. Introduction 

Systems that are driven out of equilibrium by an external field, such as temperature 
gradient or electric field, have been studied extensively in recent years. In the absence of 
a general theory, insight into their properties can be acquired by investigating simplified 
models. Studies of numerous driven models involving some conserved quantity have 
shown that their steady state typically exhibits algebraic decay of correlations [l[-|5| and 
in some cases long-range order and symmetry breaking in one dimension |6|4lO||. One 
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particular model which has drawn much attention recently is the ABC model 9|, 
This is a prototypical model for phase separation in one dimension. It consists of a 
periodic lattice of length L where each site is occupied by one of the three species 
of particles, labeled A, B and C. The model evolves by random sequential updates 
whereby particles on nearest neighbour sites are exchanged with rates 

AB^BA BchcB CA^AC. (1) 

111 

While for q = 1 the model relaxes to an equilibrium state with homogeneously 
distributed particles, it exhibits phase separation for any finite value of g 7^ 1 in the 
limit of L — 7- 00. Generically, for arbitrary choice of the number of particles of the three 
species, A^^, Nb and Nc, the model does not obey detailed balance and it relaxes to a 
nonequilibrium steady-state. A unique feature of the ABC model is that in the special 
case of Na = Nb = Nc the dynamics obeys detailed balance with respect to an effective 
Hamiltonian with long-range interactions for arbitrary value of q. This Hamiltonian 
provides a rare opportunity to gain insight into the mechanism behind phase separation 
in one dimension induced by a drive in the bulk. 

The ABC model has been considered in the weak asymmetry limit where q 
approaches 1 in the thermodynamic limit as g = exp(— /3/L) [11]. When the rate 
of approach is faster than a critical value, namely for /3 < /3c, the model reaches a 
homogenous phase in the limit of L — )■ 00. For /3 > /3c the model reaches an ordered 
phase with three macroscopic domains, each predominantly occupied by one of the 
species. The phase transition has been studied in the hydrodynamic limit by analyzing 
the linear response of the homogenous phase to small inhomogeneous perturbations. 
In the equal densities case, the transition was found to be continuous, taking place at 
/3 = /3c = 27r-\/3- The transition remains continuous for small enough deviation from the 
equal-densities condition, and becomes first order beyond a tricritical point at larger 
deviations. Since the analysis was based on linear stability of the homogenous phase, 
the full phase diagram of the model and the nature of the first order transition could 
not be explored. This would require the knowledge of the density profiles of the three 
species in the ordered phase. 

In the present paper we derive an exact expression for the steady-state density 
profiles of the ABC model for arbitrary values of average densities and /3 by solving 
the hydrodynamic equations corresponding to the evolution of the model. We use our 
results to investigate its phase diagram and the nature of the first order phase transition. 



Beyond the tricritical point we find a range of temperatures wliere botli tlie liomogeneous 
and ordered pliases are locally stable. The phase to which the model eventually relaxes 
could in principle be determined by minimizing the large deviation function. Since this 
function is known only in the limit of weak drive (/3 <^ 1) [llj and for small deviations 



from the homogenous phase 12J, we can only draw the stability limits of each phase. 

These limits define the region of parameter space where both phases are locally stable. 

The ABC model has recently been generalized to include particle-nonconserving 



processes and its phase diagram has been analyzed in the equal densities case J13L Il4| . 
The phase diagrams of the canonical (particle-conserving) and grand canonical (particle- 
nonconserving) ensembles have been shown to be inequivalent. This is in accordance 
with what is generally expected in equilibrium systems with long-range interactions. 
The study presented in this paper of the nature of the phase-separated state can be 
generalized to the case of the nonconserving ABC model with arbitrary densities. This 
would enable one to explore phenomena such as inequivalence of ensembles in a genuinely 
driven model which does not obey detailed balance [15|. 

The paper is organized as follows. We first present a brief review of the ABC model 
and previous studies of its phase diagram in section O We derive the steady-state of the 
hydrodynamic equations of the model in section [3l and express it explicitly in terms of 
elliptic integrals in Appendix A In section H] we study the resulting phase diagram of the 



model and compare it with results from Monte Carlo simulations. The low temperature 



(strong drive) behaviour of the solution is derived in Appendix B 



2. Phase diagram derived from stability analysis 

In this section we present the ABC model and review its properties and phase diagram, 
obtained in previous studies using stability analysis of the homogeneous phase. 

In studying the ordered phase of the ABC model one notes that for g < 1 the ordered 
phase is such that the domains are arranged clockwise as AA . . . ABB . . . BCC . . .C, 
and counterclockwise for g > 1. Throughout this paper we consider g < 1. The case of 
g > 1 is obtained by permutating for instance the labels of the B and C in a system 
where the drive is given by g' = 1/g < 1. 

As a result of the dynamical asymmetry, the model generically reaches a 
nonequilibrium steady state with nonvanishing currents of particles. The current of, 
say, the A particles is proportional to the rate at which they perform a full clockwise 
trip minus the rate of the counter-clockwise trip, yielding 

j^r^qNs_qNa_ (2) 

The other currents are obtained by cyclic permutation of the A, B and C labels. While 
these currents vanish in the thermodynamic limit for arbitrary average densities, in the 
special case where A^^ = Nb = Nc = -^/3, the currents also vanish for finite systems 
with arbitrary length. In this case the dynamics obeys detailed balance with respect to 



an effective long-range Hamiltonian given by 

•H (C) = 5^ 5^ -^ {AB,+k + B,Q+k + QA+fc) , (3) 

j=i fc=i 

where ^ = {C«}j=i denotes a niicrostate of the system such that Q = A,B or C. The 
operators in the Hamiltonian are defined as 

A=h «T^. (4) 

' I else ^ ' 

and similarly for Bi and Cj. The probability of a microscopic configuration is given by 
P [C,) oc q^^''\ The Hamiltonian yields a super-extensive energy which scales as E r^ L'^ 
with the system size, typical of systems with long-range interactions. 

As mentioned in the introduction, the ABC model is often considered in the limit 
of weak asymmetry , q = exp (— /3/L), where /3 is regarded as the inverse temperature of 
the model 



ll| . This rescaling of the drive with L corresponds to the Kac prescription 



for the rescaling of the temperature in long-range interacting systems 16] . It amounts 
to an effective rescaling of the energy so it becomes linear with the system size, thus 
comparable to the entropy, S ^ L. Study of this limit for the equal densities case 
revealed a second order phase transition at /3 = 27rv3 from the homogeneous state, 
where entropy dominates, to the ordered state which is dominated by the energy term. 
The ABC model has also been studied on an interval, where zero flux boundary 



condition is considered J17H19|. In that case the model obeys detailed balance also 
for nonequal densities and its steady state can be obtained using the same effective 
Hamiltonian ([3]). The steady-state density profiles of the three species in the phase 
separated state of this model has been evaluated for arbitrary values of average densities 



17|. In the special case of equal densities, the steady state of the model on an interval 
and that on a ring are related by a trivial mapping, allowing us to use the studies of 
the model on the interval as a point of reference for the present work. 

In the case of equal densities on a ring or arbitrary densities on an interval, where 
an effective Hamiltonian can be defined, it has been demonstrated that due to the weak 
anisotropy limit local density correlations vanish for L — >■ oo. Namely, 

(x,r,+i) = (X,) (r,+i) + 0(1/1). (5) 

where X, Y denote either A,B or C and ( ) denotes an ensemble averages over the 
steady-state distribution. It has been argued that this lack of local correlation is valid 
also for nonequal densities on a ring [ITJ. As a result of ([5]) the hydrodynamic equations 



171 . I20I . l2l| corresponding to this model are given by 



dp Op 

^ = /3^[p.(p.+i-p.+2)] + ^, (6) 

where r is the macroscopic time-scale and Pa{x) is the coarsed-grained density profile 
of particles of type a for a; G [0, 1]. The index a denotes the species and runs cyclicly 
over A , B and C. The conservation of particles implies that j^ dxpa{x) = N^/L = r^, 



where r^ is the average density of particle a. Since r^i + r^ + r^ = 1 it is convenient to 

express the densities in terms of two independent variables as 

111 27r Att 

(ta, tb, ^f^) = (3 ' 3 ' 3) + 2^(sin 0, sin (0 + y), sin (0 + —)), (7) 

where A^ = | Yla=A b c(^a ~ |)^ ^^ ^ measure for the deviation from equal densities and 
(/) is a phase variable. 

It is easy to see that the homogenous profile, Pa{x) = Tq,, is a solution of (|6]). Its 
stability with respect to small anisotropic perturbations revealed a critical line given by 

VI - 36A2 ^ ' 

The homogenous phase was found to be unstable at temperatures (T = 1//3) below 
this line |[ll|. Probing the region just below the critical line, infinitesimal perturbations 
around the homogenous phase were found to be stable only when r^ obey 

S (r^, rs, re) = {r\ + r| + r^) - 2 {r\ + r| + r^) < 0. (9) 

For these values the model undergoes a continuous second order transition at ([8]), 
whereas for other values of Tq, the transition becomes first order. The tricritical line, 
where the order of the transition changes, is given by S (r^, r^, 1 — r^ — r^) =0 and in 
terms of A, by 

lO8sin(30)A^-54A^ + l = 0. (10) 

The resulting phase diagram is shown in figure [1] for the case of two nonequal 
densities, defined by taking = 77r/6 as 

rA = rB = l/3-A, re = 1/3 + 2A. (11) 

The critical line and tricritical point in the figure are based on the work of Clincy et al., 
while the upper stability line is drawn based on the results presented in the two following 
sections. Note that the phase diagram is not symmetric around A = 0. At one end of 
the phase diagram, for A = 1/3, we obtain re = 1 and hence no dynamics, whereas for 
A = —1/6 we obtain the weakly asymmetric exclusion process with rA = rB = 1/2 jjj. 
While the critical line and the tricritical point can be found by expanding ([6]) near 
the homogenous solution, studying the first order transition and the stability limit of 
the phase separated state requires the knowledge of exact density profiles. In this paper 
we calculate the steady-state density profiles of the model and use them to analyze its 
complete phase diagram. An exact solution of ([6]) on an interval has been derived by 
Ayyer et al. [l7|. Following a similar derivation, we generalize their solution to the 
nonequal-densities regime of the periodic model. This allows us to study the nature of 
the first order transition phase predicted by Clincy et al. 
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Figure 1. (T = 1//3, A) phase diagram of the ABC model for two nonequal densities. 
The sohd line denotes the second order phase transition which turns into first order at 
the tricritical point (TCP). The first order transition is depicted by two stability line 
(dashed lines) where the two phases coexist. 



3. Density profiles for nonequal densities 

In this section we study the steady-state solutions of the hydrodynamic equations of 
the ABC model ([6]). Assuming time-independent solutions, we omit the time derivative 
from (E]) and integrate over x to yield 

dpa 



dx 



-/3 [Pa {Pa+1 - Pq+2)] - Ja 



(12) 



where the constants of integration, J^, are interpreted as the steady-state currents of 
particles. They can be shown to obey Ja + Jb + Jc = ^- In order to solve ( IT2l) we 
apply several transformations which are similar to those used by Ayyer et al. for the 
ABC model on an interval [l'ii|. Multiplying (fT2l) by Pa+iPa+2 and summing the three 
resulting equations yields 



^ ipAPBPc) 



Ja PbPc - Jb PcPa - Jc PaPb 
J A dpc Jc dp A 



(13) 



P dx (5 dx ' 
Integrating this equation over x yields a simple relation between the density profiles, 

PaPbPc = K -QaPc + QcPa, (14) 

where Qa = Ja/P and K is a constant of integration. One can check that this equation 
is indeed invariant under cyclic permutations of A, B and C up to a change in the 
constant of integration, K. Equation (IT^ is a generalization of the relation obtained 
for the equal densities case, where Pa{x)pb{x)pc{x) has been shown to be constant in 



space 171, |2l|. Using fll4p in conjunction with Pb = ^~ Pa~ Pc allows us to express p^ 



in terms of pc as 

- (Pc - PC + Qc) ± \l{pl -PC + Qcf + ^QapI - ^Kpc ^ ^ 

PA = ^- (15) 

^Pc 

Inserting this expression back in (fT2l) for a = C yields an explicit equation for pc, 
dpc 



Q^ ±\l{pl~Pc + Qcr + ^QAPl-^Kpc. (16) 

The plus and minus signs correspond to the two halves of the ring around the maximum 
of pc{x). Taking the square of this equation and writing it in the rescaled variables 

t = 2/3 X and y{t) = pc{x) we obtain 

V(t)2 + [/x,Q,,Q,(l/(t))=0, (17) 



2' 



where 



1^2,, ^2 , 2K + Qc 2Qa + Qc 2 Ql 



Uk,q.,qc (y) = "8^(1 - yy + — ^y - ""^ ^ y' - ^. (is) 

Equation (fT7|) can be viewed as an equation of motion of a zero-energy particle with 
mass 1 in a quartic potential. Equation (TTTll and the derivation below can be written 
in terms of either of the three species by cyclic permutation of A, B and C . In the 
nonequal densities case the quartic potential changes under this permutation, yielding 
a different profile for each species. 

Depending on the values of K, Qa and Qc the potential may have two, three or four 
real roots, depicted in the (i),(ii) and (iii) lines in figure [21 respectively. The four roots 
of the potential, denoted as {a, b, c, d} can be shown to obey 0<a<b<c<l<d. 
In this case the particle oscillates between b and c which is the only physical trajectory. 
This is because we require that both < y{t) < 1 and U{y) < 0. The case of three roots, 
when b = c, yields a constant trajectory in time which corresponds to the homogenous 
solution, pc{x) = y{2(ix) = re- The case where there are only two real roots does 
not correspond to any physical solution. The manifold which defines the region of 
{K, Qai Qc}-space where the physical solution resides is thus obtained by inserting the 
homogenous solution, Pa{x) = r^, into (fT2|) and f lT^ as 

Qa,h = ra ira+2 - ra+i) , Kh = TATBrc + QArc- Qc rA- (19) 

The trajectory of ( TTTj) between b and c for b < c yields the ordered profile of the 
ABC model corresponding to given values of K, Qa and Qc- In order to relate these 
parameters to the original parameters of the problem, /3 and r^, we examine the period 
of oscillation of the particle between b and c, which we denote as O. The periodic 
boundary condition of the profile, pc{x + 1) = Pc{^)i imposes a constraint on the 
solution of the form, ymii + 2/3) = ymif) or equivalently = 2/3 /m. Here, the positive 
integer parameter m corresponds to the number of times the particle oscillates between 
b and c in a time interval of length 2/3. We argue in the next section that only the 
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Figure 2. The effective potential, C/(?/), for ta = tb ^ re and some values of K, Qa 
[Qc = 0). (i), (ii) and (iii) represent the case where potential has 2,3 and 4 real roots. 
They correspond to the case of the model has no physical solution, a homogenous 
solution and an ordered solution, respectively. The four roots of the potential in (iii) 
are denoted on the graph by {a, &, c, d}. (iv) depicts the limit of if ^ 1, Q <C 1 which 
corresponds to the low temperature limit (/3 ^ 1) . 



m = 1 solution describes the ordered steady state of the model. The periodic boundary 
condition may be written as 



dx 



^^ di_ 

2/3 



m 



dy 



'0 ^0 ^/^ 1^ Jb \/-'^UK,QA,Qciy) 

An additional constraint on ym{t) comes from the total number of C particles, 



re 



dxpc{x) 



'^ dt 



ym{t) 



m 



,0 2/3--^^^ (5 h ^-2UK,Q,,Qa{yy 



ydy 



(20) 



(21) 



The third constraint is obtained by dividing (fT2l) by pa and integrating the result over 
X using periodic boundary conditions. For a = C the result yields the condition 



rs -ta 



dxp 



^^ dt 



^C Jo Jo ^P P Jb \/-2UK,QA,Qc\y) 

which is related to the difference between r^ and r^ and the consequent current of C 



y ^dy 



particles. In Appendix A we provide an analytic expression for (I20|) - (l22|) using elliptic 
integrals. 

In the following section we will analyze the phase diagram which arises from the 
solution above. For simplicity we restrict ourselves to the two nonequal densities case, 
fA = ''^B^fCy which yields the same qualitative behaviour as the more general three 
nonequal densities case. For r^ = Tb y^ re we find that Qc = Jcl P = 0. This simplifies 
the form of the effective potential 0181) and leaves us with only two constraints, (I2U]) and 
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Figure 3. The ordered density profile obtained from (P5)) and (IT5|) for va — tb ^ re 
with A = 0.1, 13 = 25 and to = 1. 



( pTj) . which take the form of 

Q = 2l3/m = 2xK{l/k)/k, 

1 + a_/a+ 



1 , 1 

re = - + 2A = — 
6 a;_ 



(5k/mx 



Y{{al/alk'',l/k) 



(23) 
(24) 



Here a±, x and k are functions of K, Qa given in Appendix A and K{k), Il{n, k) denote 
the complete elhptic integral of the first and third kind, respectively, whose definition 



is found in [Appendix A| as well. The profile of the C particles is expressed by inverting 



the equation x = Ldt'=^ P , ^ 

^ Jo 2/3 Jb ^-2U{y') 



as 



Pc{x) 



1 + sn (2/3x/x, k) 



Ua 



(25) 



a_sn (2/3x/x, k) 

where sn is the Jacobi's elliptic function [22|. The dependence of the profile on m is 
hidden in the value of K,Qa which set k,x,a±. The resulting m = 1 profile for a 
specific value of /3 is shown in figure [31 In the section below we study the behaviour of 
this solution and the resulting phase diagram. We also examine its behaviour at low 



temperature (1 ^ /3 ^ L) in Appendix B and find that it conforms with our physical 
understanding of the model. 

It is interesting to note that the hydrodynamic equations of the ABC model ([6]) 
can be solved by considering a moving steady-state solutions of the form p^ {x, r) = 
Pa{x + VT). However, such solutions did not appear in the numerical relaxation of 
([6]) as well as in Monte Carlo simulations. This may mean, although remains to be 
proven, that moving solutions are unstable stationary solutions of the ABC dynamics. 
We therefore restricted our derivation to case oi v = 0. 



4. Phase diagram for nonequal densities 

In this section we examine the behaviour of K, Qa as we change the values of /3 and r^ 
for the two nonequal densities case. Their values are obtained by inverting the integral 
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conditions given in fl20p - (l^T]) . which are written in an exphcit form in fl23p - (12^ . 

For all values of /3 the hydrodynamic equations Q have a stationary solution of 
the form Pa{x) = r^, which corresponds to the homogenous values oi K = Kh and 
Qa = Qh,a given in (IT9|) . As discussed in section [21 this solution becomes unstable 
below the critical line, 13 > (3c = 27rv3/A/l — SGA^. In this regime we expect to find an 
ordered solution. Figure H] displays T = 1//3 computed according to (123|) where Qa is set 
for a given value of K through (12^ . For a small values of A, in figure H^, we find a second 
order transition at T = 1/ fie between the homogenous phase and the m = 1 ordered 
phase where K < K^. This behaviour persists up to the tricritical point ([9]), which in 
the two nonequal densities case takes the simpler form of A = 1/ (3 + 3v^) — 0.122. In 
figure IDd, we see that beyond the tricritical point, the m = 1 ordered phase appears also 
at T > l//3c- This is because the relation between K and /3 under a fixed value of A 
is non-monotonic. As a result the model is expected to undergo a first order transition 
between the two phases at value of /3 between the two stability limits. The discontinuity 
in K at the transition implies that this is a transition from a homogenous state to an 
ordered state with a finite amplitude of modulation. 

In order to compute the first order transition point one has to know the full large 
deviation function (LDF) of the ABC model, which is not known. We may still draw the 
stability limits of the two phases defined by the critical temperature and the minimum 
of K{P) in the ordered phase. The resulting stability lines are shown in figure [H for the 
case of two nonequal densities. 

In figure [5] we examine the first order transition using Monte Carlo simulations. 
The algorithm for the simulation is straightforward. At each step a site is selected at 
random and an exchange step is attempted where the particle in the chosen site may be 
exchanged with its neighbour to the right with probability given by ([1]). We measured 
the parameter (paPbPc) aiid compared it to that obtained from the hydrodynamic 
solution by integrating (HM over x. In the simulation it was measured by counting the 
number of ABC triplets in the lattice after each L exchange attempts and averaging 
the result over many such time steps. 

In figure [5] we plot the simulation results for different values of temperatures around 
the first order transition and for various system lengths. For each value of T the 
simulation is initiated in the fully ordered phase and run for a time period which was 
sufficient to observe transitions between the two phases. The number of ABC triplets 
is averaged over the entire second half of the simulation where the system is unaffected 
by its initial state. The figure displays a first order phase transition, smoothen by 
finite size effects. The transition occurs below the critical point as suggested by our 
analysis. Near the transition point we observe slow fluctuations of the system between 
the two phases, as depicted in flgure [6l This implies that the flgure [5] might contain 
some errors near the transition point due to insufficient sampling time. We do not 
expect the transition point in the L — )■ oo limit to obey Maxwell's construction since 
the horizontal axis is not the conjugate variable of T. The latter can only be derived 
from full LDF of the model. Figure [5] also displays a good agreement with the theoretical 



11 



K/K, 



K/K, 




Figure 4. The temperature, T = 1//3 = 2/mQ{K,QA), computed as a function of 
K for m = 1,2,3 for two nonequal densities, r^ = rB ^ tq- Qa{K,A) is set by 
([24]) . The dashed hue denotes the critical point (3 — 2n^/3/\/l — 36A^. (a) and (b) 
are calculated for A = 0.12 and A — 0.14, respectively, depicting the case where A is 
below and above the tricritical point. 



values for (paPbPc) above and below the transition point, which confirms the validity 
of the mean- field approximation ([5]) . 

Figure [6] depicts the fiuctuations of the system between the ordered and disordered 
phases for L = 4800 at a temperature close the first order transition point. The figure 
shows significant and long-lived fiuctuations around the ordered phase. A thorough 
investigation of their nature showed that they are not related to any known meta-stable 
state of the model and that they decay as the size of the system is increased. We avoid, 
however, using larger systems since they would require much longer simulation time to 
display transitions between the two phases. 

In figure m we find that for l/2/3c < T < 1/Pc the m = 1 is the only stable solution, 
whereas the m > 1 profiles are unstable. At lower temperatures the latter become 
stationary states and may theoretically be the ground state of the model. In the equal- 
densities case, this possibility has been ruled out by showing that the m = 1 profile 
has the lowest free energy for all T < T^ |l7|. For nonequal densities, a similar analysis 
would require the knowledge of the full LDF of the ABC model. Here, however, the fact 
that the m = 1 is the ground state of the model can be argued by noting that the m*'^ 
solution corresponds to an ordered state with particles segregated into 3m domains. 
Since lower temperatures (stronger drive) favour segregation, it is natural to assume 
that the most segregated state, m = 1, remains stable for T < l/2/3c. We therefore 
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Figure 5. The density of ABC triplets as measured in simulation as a function of 
T = 1//3 in comparison with the theoretical values (solid line). The simulation was 
performed for the two nonequal densities case with A = 0.14 > A^cp — 0.122 and 
various values of L. The dashed line denotes the value of /3 = 19 for which the time 
evolution of {pApBpc) is plotted in figure [6] 
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Figure 6. The density of ABC triplets as a function of time given by the number of 
Monte Carlo sweeps for L — 4800, A = 0.14 and /3 = 19. The upper and lower dashed 
lines denote the theoretical values of the homogenous and ordered phases, respectively. 
The dotted line denotes the unstable ordered solution which, as expected, does not 
appear in the simulation. 



consider it to be the steady-state solution of the model throughout the ordered phase. 
This argument is supported by Monte Carlo simulations of the model, where only the 
m = 1 profile was observed below the transition point. 



13 
5. Conclusions 

In this paper we liave derived an exact expression for steady-state density profile of 
the ABC model on a ring. The derivation is based on a hydrodynaniic description of 
the model which has been suggested to be exact in the thermodynamic limit [IT]. The 
solution is valid for the case where the average densities of each species are not equal 
and the model is thus out of equilibrium. Using this solution we have studied the first 
order phase transition whose existence was suggested by Clincy et al. The transition is 
of first order since there is a finite range of temperatures where the model has two stable 
phases. The transition point is located where the large deviation function of the two 
phases is equal. However, since this function has not yet been found we can only draw 
the stability limits of the two phases which define the coexistence region. Monte Carlo 
simulation of a specific point in parameter space confirmed that the first ordered phase 
transition occurs within the coexistence region, above the previously known critical 
temperature. The simulations also yielded a good agreement with the hydrodynaniic 
solution in regions where only one phase is stable. 

The present study opens the door for future studies of the ordered phase in the 
ABC model with arbitrary values of average densities. We plan to employ the solution 
obtained here in order to study the inequivalence of ensembles in the nonconserving 
ABC model with nonequal densities |15| and compare the results with those previously 



obtained for the equal-densities case J13l . Il4 
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Appendix A. Analytic expression of the mean-field solution 

In this appendix we obtain an analytic expression for the relation between the 
parameters of the effective potential K, Q^, Qc and the parameters of the model, /3, r^. 
We derive the explicit expression for the general case and then present its simpler form 
for the two nonequal densities case, r^ = Tb ^ re, mostly used in this study. 

Appendix A.l. Three nonequal densities 

In section |3] we mapped the mean-field dynamics of the ABC model to the motion of a 
particle in a quartic potential, 

TT I \ '^ 2/1 \2 , 2i^ + Qc 2Qa + Qc 2 Qc /A 1^ 

UK,QA,QcKy) = -^y K^-y) + — ^ — y ^ — y --^- (a.i) 

The parameters of the potential K, Qa, Qc are linked to the parameters of the model, P 
and Ta, through three conditions (l20l)-(l22l) which involve integration over the trajectory 
of the particle. 
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We begin with fl20p which can be expressed in terms of the incomplete elhptic 
integral of the first kind defined here in the Jacobi form |2^ 



J ^{l-z^){l-k^z'')' 



Following a similar derivation as in [17| we introduce a Mobius transformation that 
maps the roots of the potential, 

UK,QA,Qciy)=^ "^y ^{a,b,c,d}, (A.3) 

onto the poles of the integrand in (1A.2P , 

{a,6,c,d}^ {-1,-1/A;, l/A;,!}. (A.4) 

The transformation is given by 

.= f(y) = ^Z^ ^ + ^-^ ^Ay + 5 



where 



and 



±(bc — ad) + ^y(a — b)(a — c)(b — d)(c— d) 

"± = h f , ^\ — i7rT~\ ' ^^-^^ 

oc[a + d) — ad[o + c) 



L-t m- yu ^ab~2ad+cdj 

1 ry (b n i^-^)^^-'^) \ 
l-a+ \0-aab-2ad+cd) 



(A.7) 

The parameters a_,a+, A, B, C, D and k are functions of K^ Qa, Qc through a, b, c and 
d. Let t{y) be the time it takes the particle to move from b to y. Using the transformation 
above it may be expressed as 

y f{y) 

f dy' f dz 

t=2 . ^ = X ' 



-2U{y') J V(l-22)(i_p^2) 

" — i/fe 

= x[Fil/k,k) + F{f{y),k)], (A.8) 

where 



2^{C'~A^)iC^-k^A^) 
"^ AD-BC ^ ' 

Equation (l2Up . whereby the period of oscillation obeys = 2/3 /m, can be expressed 
by setting f{y) = 1/k in (lA.Sp . For that end we first notice that the integral in this 
equation can be brought to a simpler form in the new coordinates w = kz, 

i/fc 

F{l/k,k)= I , "^^ = 

J J(l-z^)(l-k^z^) 
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where K (l/k) is the complete elhptic integral of the first kind. Using this form the 
condition of = 2/3/m can be written as 
2mx. 



'-K{l/k) = f3. 



k 
The second condition ( 12T]) can be written in the form of 

"2^ ^ /-C ydy 



(A.ll) 



1 f^^ m r 



mx /"fe —B + D2; 



-2f/(z/) 
dz 



/3 J^^_ A-Cz ^{l-z^){l~k^z^) 

mn /"^ 2 (aD-BC 

~TJo AC 



(A.12) 



-DA 



(i^; 



2mxD 

/3C 



1 






v/(l-22)(l-fc2^2) 



^'P* 



Here n(?T,, x, rra) is the incomplete elliptic integral of the third kind 

X 

Il{n,x,k) 



(1 - nz^) ^{l-z^){l-k^z'^y 



22l | defined as 

(A.13) 



The third integral condition ( l22l) is given by replacing A with B and C with D in ( 1A.12I) . 



2/3 






y\ly__ 



■mx: 



-A + Cz 



dz 



13 J_i_ B-Dz ^(l-z2)(i_p^2) 

2??2xC 
/3D 



(A.14) 



AD\ fD^ 1 , 
1-^r^ in — ,-,A; 



^•P^ 



As in ( lA.lOp , these results can be written in terms of a complete elliptic integral using 
the transformation w = kz. 



i/fc 



U{n,l/k,k) 



dz 



nz^)^ 



k^z^) 



dw/k 



(l-1^)\/(l-f^)(l 



w^ 



k ^k^'k' 



(A.15) 



Here H(ra, k) is the complete elliptic integral of the third kind. The two integral 
conditions (12T]) . (122|) are thus given by 



D 

C 

15 



\^ l-BC/AD ^ 

(ik/mx 



C 



k'^A^' k 



I- AD/BC ( D'^ 1 



(5k/mx 



k'^B^' k 



re, 



Q 



c 



(A.16) 
(A.17) 
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In order to obtain K{P, Tq,), Qa{P, fa) and QcW, ^a) one has to invert (lA.lip . (lA.16p 
and flA.17p . Once this is done the profile can computed by inverting (lA.Sp . The result 
may be expressed using the Jacobi elliptic function, sn{x,k), defined by the equation 
F (sn {x, k) ,k) = x. The profile of C particles is then given up to translations of x by 

Note that the dependence on m is hidden in k, x, a±. The two other profiles, Pa{x) and 
Pb(x), are obtained from ( !T5|) and Pa{x) + pb{x) = 1 — pc{,x). 

Appendix A. 2. Two nonequal densities 

For convenience we write explicitly the solution for the two nonequal densities case, 
which is studied extensively in this paper. The form of the solution in this case is very 
similar to that obtained for equal densities in [17i |. 

When ta = tb ^ re it easy to see from symmetry that Qc = Jc/ P = 0. As a 
result, our solution depends only on two parameters K, Qa which are set by /3, A. The 
latter is defined by 

rA = rB = 1/3 - A, re = 1/3 + 2A. (A.19) 

The effective quartic potential which governs the motion of the particle is of the form 

^K,Q.,o {y) = -ly'ii - yf + f y - ^y' (A.20) 

and its four roots are thus {0,b,c,d}. The Mobius transformation from these roots to 
the poles of the elliptic integral are obtained by setting a = in the form presented in 
the previous section. More explicitly it is given by 

z = f{y) = ^^^, (A.21) 

a^y + 1 



where 



and 



±bc + ^hc {d -c){d- h) 
a± = ^ , (A.22) 



* = i±^. (A^23) 

1 — a+0 



Note that this form is identical to the one defined in the equal densities case [Ij 

We now express using elliptic integrals the two conditions that determine mapping 
between K, Qa and /3, A. The first condition is identical to (lA.lip . 



2mKK (1/fc) /k = /3. (A.24) 

where here x takes the simpler form of 

2(0:4- + a_) 



X 



a/(1 — a;+6)(l — q;_|_c)(1 — a^d) 



(A.25) 



The second condition (121 p is given in this case as 
1 



a. 



'^ l + a^/a+ ^ ( a^_ 
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re = ^ + 2A. (A.26) 



fikj Km ya^fc^' k 

The functions ir(/3, A), (5a(/3, A) are obtained by inverting ([Oil) and ([06D . As 
in the previous section, we use the resuh to express the profile of the C particles, 

= 1 + sn (2/3x/x, fc) 2 

«+ — a_sn (2/3x/x, fc) 

The two other profiles, Pa[p^) and Pb{x), are again obtained from (fT5|) and pA(a^) + 
Pb(x) = 1 -pc(a;). 

Appendix B. Asymptotic behaviour at low temperatures 

In this section we study the behaviour of the hydrodynamic solution of model at low 
temperature (L — )■ oo, 1 ^ /3 ^ L) for the case of two nonequal densities, va = tb y^ re- 
This form will be especially useful in future studies of a generalized ABC model with 



nonconserving dynamics [15|. The T = limit corresponds to i^ = Qa = 0. Inserting 
this into (1X221) . (1X231) and (lAl25l) we obtain 

k = l, a± = 1, x = 4. (B.l) 

For k = 1 the elliptic integral in (1A.24P diverges, corresponding to the limit of /3 — t- oo. 
The elliptic integral in ( 1A.26I) diverges as well, maintaining a finite value of re- 

To study the low temperatures behaviour we assume that i^(/3. A) and (5a(/3, A) 
vanish for exponentially to leading order as 

K^e-^^\ QA^e-^^\ (B.2) 

In order to find the coefficients 71 and 72 we first examine the behaviour of (lA.24p . 

Here we analyze only the behaviour of the m = 1 solution, which is considered to 
be the ground state of the model (see section H]). We expand x and k around their value 
at T = as 

x = 4 + >€i{K,Qa) k=l + ki{K,QA) (B.3) 

where Xi and ki are functions whose form is not written explicitly in order to avoid 
lengthy expressions. They obey < Xi ^ 1 and < A;i ^ 1. Expanding (1A.24I) to 
leading order in these functions J23| we find that 

ln(A;i) = -/3/4 + 0(l) (B.4) 

The condition coming from the average density (1A.26I) involves an additional function, 

a^/a+ = l-ai{K,QA) (B.5) 

which can be shown to obey fci ^ ai ^ 1. Expanding (IA.26P in ai and fci to leading 



order |23| we obtain that 

ln(A;i) - 3 In(ai) = 12A In(ai) + 0(1). (B.6) 
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Solving (IB.4p and (1B.6P to lowest order in K and Qa yields 

^-^ ' ^A ~ I ^_^(i/3+2A) A<0 ' ^^^ 

which agrees with form of i^ ~ e~^/'^ found in the equal densities case [19]. 
In order to interpret this result we observe the behaviour of ( TT^ whereby 

Pa{x)pb{x)pc{x) = K -Qa Pcix), (B.8) 

and hence 

dxpA{x)pB{x)pc{x) ~ e-^'"'"(^^'^+2^) = q^MNA,N^'^c)_ (B.9) 



'0 

In the limit of /3 — t- oo the number of triplets of ABC is governed by the probability 
of an event where in the fully separated state two particles of different species meet in 
the domain of the third species. This probability scales as q"^^^(^A,'^B,Nc)^ because this 
event occurs with the highest probability in the smallest domain. 
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